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We discuss the implications of the eikonal amplitude on the pair production 
probability in ultrarelativistic heavy- ion transits. In this context the Weizsacker- 
Williams method is shown to be exact in the ultrarelativistic limit, irrespective 
of the produced particles' mass. A new equivalent single-photon distribution is 
derived which correctly accounts for the Coulomb distortions. As an immediate 
application, consequences for unitarity violation in photo-dissociation processes in 
peripheral heavy-ion encounters are discussed. 

1 Introduction 

The S matrix describing electron scattering at ultrarelativistic pointlike charges was shown 
to be determined by the gauge phase leading to the Dirac equation represented in the 
temporal gauge |I] . We found that it naturally exhibits the same form as the well known 
eikonal expression, as is expected by Lorentz invariance. 
The gauge phase leading to the temporal gauge reads 



The consideration of asymptotic states corresponds to sending the upper bound of the 
integral to infinity. The infinite time integral over the scalar part of the electromagnetic 
potential in the exponent has to be understood as the principal value of the integral. It 
can be decomposed into a finite term and an infinite quantity, expressing the familiar 
divergence of phases in Coulomb scattering. The infinite term can be removed by a gauge 
transformation which converts the Coulomb boundary conditions of the original problem 
into a modified short-range potential allowing for asymptotic plane wave solutions. 
The ultrarelativistic limit of the gauge transformed potential A' (x) reads |2|, |, [|, [1J 



and hence we obtain for the S operator in coordinate space (7_ =70 — 73 is the Dirac 
matrix structure of the interaction) 
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The obtained S operator is a unitary operatorb due to its conformity to the gauge phase. 
It agrees with the first term of the Magnus expansion of the time-evolution-operator 
f|, since the considered gauge-transformed interaction was assumed to be compressed to 
infinitely short times. 

This result was proven to be of completely perturbative nature [1]]. Note, however, 
that the perturbative derivation did not require the explicit Fourier transform of the 
transverse part of flU), which is an ill defined object. For that reason the deduction of the 
small-coupling limit (Za — > 0) of @ in momentum space must be treated with special 
care. 

In a rigorous distributional sense it can be defined as the weak limit A — > || 

dVe-^- lnxi = lim4vr - n6 2 (k ± ) In (j^)) (3) 

with A = 07/7, fi = 2/e c . 'Weak limit' means that the limit A — * has to be taken 
after having integrated the expression with a test function. The second term on the RHS 
arises from the gauge transformation applied to the potential and thus accounts for the 
Coulomb distortions. 

Discarding the second term and taking the limit A —>■ directly would accidentally 
yield an expression for the Fourier transform of (|IJ), being identical to the Fourier trans- 
form of the ungauged potential in the limit 7 — > 00 
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/ d 4 x e lkx A = -{2n) 2 Za5{k_)-j 

J K 1 



This error is made if one intends to extract the correct small-coupling limit from a 
naive Taylor expansion of the Fourier transformed T matrix whose linear term reads 
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\im n T(k) « (27r) 2 5(k^)iZa T ^u(p)^f-u(p) (4) 



Here u(p) is the electron unit spinor, p and p' are the initial and final momenta of the 
electron, repectively, k = p' — p. 

Since, however, Taylor expansion and Fourier transformation do not commute in this 
case, the Taylor expansion of the Fourier transformed T matrix for this purpose is not 
justified.^] The correct small-coupling limit of the scattering amplitude in momentum 
space can thus not be found by a naive Taylor expansion of the Fourier transformed T 
matrix and does not agree with first-order perturbation theory. According to (|f) this 
is simply based on the fact, that the gauge transformed potential correctly accounts for 
Coulomb boundary conditions. 

In the following we want to investigate, how the correct treatment of Coulomb bound- 
ary conditions in all orders of perturbation theory influences the cross section of the 
scattering process. 

-'^Note, that with the above mentioned distribution-theoretical precautions, it is possible to obtain the 
correct result via Taylor expansion [Q. 
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2 Implications on the cross section 



We consider the exact amplitude for electron scattering at an ultrarelativistic pointlike 
charge, moving in +z direction. We state it in terms of the invariant squared momentum 
transfer t w — k\ for 7 — ► oo0 

A = 2-K5{kJ)F p ,, p {e- iZalnx ^ - l)u(p')^u{p) 

= -i8^Zatf(fc_)iS=^ e ^ a ^-'/%(^7-«(p) (5) 
t 1 [iaZ ) 

F p i p abbreviates the Fourier transform with respect to the transverse coordinates, taken 
at the momentum kj_ = (p' ± — p±). The cross section for this scattering process is found 
to be exactly the Mott formula for Coulomb scattering of ultrarelativistic electrons at a 
static source, Lorentz-transformed into the electron's rest frame. Such kind of agreement 
between the exact result and the first order perturbation theory is also found in the 
nonrelativistic case, known as one of the peculiarities of the Coulomb field. 

The well established eikonalization of the scattering amplitude and thus the reduc- 
tion to Mott's result imply, that in the high-energy limit the electron and the positron 
Coulomb scattering cross section become identical. This behaviour of the cross section at 
ultrarelativistic energies is required by the Pomeranchuk theorem ||. 

One can draw an analogy to pomeron exchange in hadron physics: The scattering 
process can be described in terms of the single exchange of an "effective photon" according 
to the modified potential 

V (x) = V 3 (x) = 5{z - t) {^) lZa ~ l) (6) 

In the field of two ultrarelativistic colliding pointlike nuclei, the exact scattering am- 
plitude was shown to retain the structure of the second-order perturbative result, due to 
the causal decoupling property Jl|. Each interaction can be described by the modified 
potential (|6|). In the following we consider the symmetric collision of two ions with charge 
Ze, the extention to asymmetric collisions is trivial. 

Accounting for both time orderings, the amplitude reads 



Apl = I ^±F k , p (e- lZaln ^ - 1)F P , fc (e- iZoln *i - l) e *^x-^)-5 



_/ n -ao. • k± + 7 m A 
p + P- — kj_ — m z + le 

,s -ol ■ (P± +P'± ~ k±) +70^ * f \\ 

u ( p ) ~ 7 1 — rv2 — ^—^- u (p) (7) 

P-P+ — (p± + p j_ — k±) z — m z + te ' 



2 Note the striking similarity between and the nonrelativistic (Schrodinger) amplitude 



2Psin^fr(l-^ 



with the squared momentum transfer being proportional to sin 9/2 
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Here the trajectory of one ion is shifted by the impact parameter b, which is accounted 
for by the factor e l( - p ' J -~ k± ^ b . 

We now use the crossing invariance of the amplitude to apply the obtained result 
to electron-positron pair production. The initial electron four momentum p has then 
to be replaced by the negative final positron momentum p — > — pp. The final electron 
momentum will be denoted by p' = p e . With (|7|) we obtain for the pair production 
probability 



da 
dR) 



— — — = \A% 



tot 1 2 



m 



d 3 p e md 3 pP 



p'p\ 



(2ir) 3 E e (2tt) 3 Ep 
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-a± ■ (~P P ± + p" e ± ~ k±) + 7 m 
—p e _pp + — (—pP L + p e ± — k±) 2 — m 2 + ie 

— t~3 — 7j2 r~ -iM-r) 

-a ± ■ {-pP L + p e ± - fc'j.) + lom 
—p e _pp + — (— + p e ± — k'±) 2 — m 2 + ie 



7-1^— pf 



The integration over the impact parameter yields the pair production cross section. 
Due to the 5 2 (k'j_ — fcjj-function occuring in the b integration, a further momentum 
integral can be perfomed, leaving 
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m 



d 3 p e md 3 pP f d 2 k 
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m 2 + ie 



7_w(-p p ) 



(9) 



Thus, upon integration over the whole impact parameter plane, the phases in the in- 
dividual scattering amplitudes (see @) cancel. Consequently, in the limit 7 — > 00 the 
cross section is found to reduce to the lowest-order perturbation theory, the two-photon 
result. This behaviour does not only naturally explain [|l(J the experimentally observed 
quadratic dependence on the target's and the projectile's charge 0, but also implies, that 
no asymmetries should occur in electron and positron spectra. 

Equation is strictly valid only for pointlike ions. The focus on electromagnetic 
reactions in peripheral heavy-ion collisions implies a restricted range of impact parameters 
with a lower bound b = ta + r_B, ta and tb being the radii of the ions. Therefore in 
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experiments which are triggered on peripheral collisions, effects of the Coulomb distortions 
described by the phase in (|5[) will be visible. 

The eikonal approximation (and thus the cross section) is known to become energy- 



independent in the ultrarelativistic limit [IT]. This dependence is restored by accounting 



for the correct transverse momentum range, which is restricted by the validity of ([!]). This 
condition reads |12|, |J 

\k±\ > - (10) 
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Such a low energy cut off is also necessary to cure the divergence in (|?p. 



3 Equivalent Photon Approximation 

We intend to study the behaviour of the cross section, both impact parameter dependent 
and impact parameter integrated, in the Weizsacker- Williams method of equivalent pho- 
tons. This approximation uses the similarity between the fields of a fast moving charge 
and a swarm of real photons moving in beam direction. It approximately correstponds to 
the first order Born approximation in the temporal gauge: Only the transverse part of the 
interaction is considered - the longitudinal part is suppressed by I/7 2 - and the vertex 
function is evaluated on-shell at k 2 = 0, i.e. for an assumed real photon. Rewriting the 
exact cross section in terms of the real photon cross section, the whole information about 
the scattering potential, which can be the retarded Coulomb potential or the modified 
potential (0), respectively, is then contained in the distribution function of the equiva- 
lent photons n(oS). Roughly speaking, this photon distribution function is determined 
by the squared absolut value of the Fourier-transformed potential (in temporal gauge). 
The obvious advantage of casting the exchange of effective photons according to (H) in 
the Weizsacker- Williams form is, that any difference between the second-order perturba- 
tive result and the exact calculation will be solely generated by differences between the 
equivalent photon distributions. 

The Weizsacker- Williams approximation is applicable if the exchanged momentum 



meets the following conditions |L3 



^<|A;_L|<m (11) 
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and 

|A;±|<^<m7 (12) 

The upper bounds mainly stem from the requirement, that \k 2 \ is negligible compared 
to 2p ■ k > m 2 , such that the intermediate propagators of the scattered particles can be 
approximated by those describing the interaction with real photons. The particle's rest 
mass in (|TT| ) is however a conservative upper bound and the equivalent photon method is 



not strictly invalid for \k ~ m . Note, that for the approximative calculations in [14 



the transverse mass of the scattered particle was taken as the upper bound for (JXTJ) . 

Since the exact amplitude takes the eikonal form, we must point out the following: 
The expansion of the ultrarelativistic scattering amplitude in powers of the transferred 
momentum yields, as the leading term, the eikonal expression (describing the minimal 



deflection from the initial straight-line trajectory) |15| , |i6| . Its perturbation-theoretical 
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derivation requires that the quadratic terms k 2 are negligible relative to the terms 2pi ■ k in 
the denominators of the propagators, where k is any partial sum of the internal momenta 
17| . The exact validity of the eikonal formula at infinite energies therefore shows, that 
the transferred momentum \k 2 \ does not exceed m 2 (irrespective of the value of m). This 
agrees with the theoretical observation, that at high energies particles are predominantly 
scattered into a cone with opening angle 9 ~ I/7, corresponding to momentum transfers 



I A- 2 



The main contributions to the cross section are thus expected from spatial 



distances larger or equal the Compton wavelength of the particle. f] 

Moreover, the longitudinal part of the interaction vanishes identically in the limit 

7 — > 00. Hence, the applicability conditions of the Weizsacker- Williams approximation 

are trivially fulfilled in the limit 7 — > 00 |J 

We have the freedom to apply this method to the interaction of the electron with both 

nuclei, giving two possibilities (see Figure [I]) 



n A (co) x Compton x n^co) Bremsstrahlung x n^co) 




a) 





b) 




Figure 1: The two possible distinct processes, that can be used to describe electron- 
positron production in heavy-ion collision. One or both ultrarelativistic ions can be re- 
placed by an equivalent photon distribution. If the bremsstrahlung process b) is calculated 
in one ion's rest frame, the electron must be assumed ultrarelativistic, to yield agreement 
with a). 

The two possible calculation schemes (Figure [I]) agree, since in b) the bremsstrahlung 
emission and the scattering at the external potential decouple. This is due to the fact, that 
the region in which the ultrarelativistic electron "feels" the external field is assumed to be 
pointlike and any frequency of the emitted photon is " soft " compared to the timescale of 
the scattering. Coulomb effects arise, if one explicitly accounts for the finite interaction 
time, either in the scattering process by correcting the eikonal formula or by keeping 
the eikonal amplitude for the scattering process but assuming a Rutherford-deflected 



trajectory for the photon emission | 13fl . Corrections to the eikonal formula account for 
higher orders of e.g. the Magnus expansion ||, which is an expansion in the interaction 
time t around the instantaneous interaction r ~ I/7 —>■ 0. In general these Coulomb 



3 From the asymmetry of electron and positron spectra produced in S(200 GeV/n)+Au collisions , the 
mean transverse distance from the target ion was deduced to be approximately two Compton wavelengths 
||. The collision energy corresponds to 7 ~ 10 in the center of speed system. 

4 Just as the eikonal formula, the Weizsacker- Williams approximation can be viewed as the leading 
term of an expansion in powers of k 2 /m 2 flg| . The validity of the eikonal expression then automatically 
implies the validity of the Weizsacker- Williams method. 
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effects vanish, if the energy of the emitted photon is too small to resolve details of the 
scattering process, and the recoil of the electron is negligible. 

To apply the Weizsacker- Williams method to the bremsstrahlung photon, the recoil of 
the bremsstrahlung photon must, however, be assumed negligible. The small momentum 
transfer is in turn ensured by the eikonalization of the scattering process. 

The equivalent single-photon distributions n^/^o;) of the ions A and B, are deter- 
mined from the effective potential (|6|). The photon distribution reads 

1 f m 

n(u) = — ^ — / k±dkx 

The lower bound of the integral is taken from the condition fllOD . The upper bound, the 
electron rest mass, is imposed by ([IT]). The prefactor arises from properly rewriting the 
cross section @ in terms of the real photon cross section (i.e. the Compton cross section) 
and photon distribution functions. 

This photon distribution coincides with the equivalent-photon distribution obtained 
from the Coulomb potential to logarithmic accuracy |T3|] and is not changed by Coulomb 
effects. 



k i irZa 




X-iZa 



T(—iaZ) 



T(iaZ) 



In . 

7VLU \ UJ 



(13) 



4 Impact parameter dependent cross section 

The impact parameter dependent equivalent photon method for the exact calculation, 
using the modified potential (Q), can be derived similarly to fll9| . To this end we have 
to modify the integrands in (|j) such that they account for the limited momentum range 
((ED and (0)). 

The cut off of low transverse momenta according to (|10|) can be achieved by the 
following replacement in (|5])[] 

* = -*!- — - k l (14) 

T 

This substitution suppresses small transverse momenta less strongly than the strict cut 
off at k± — oo / j. It assumes the sufficient accuracy of the classical eikonal amplitude for 
I/7 in the near vicinity of I/7 = 0, which is guaranteed by the possibility of continuous 
extentions of the eikonal formula towards finite 7 and large t. A Yukawa-type damping 
of transverse distances corresponding to the cut off transverse momenta yields additional 
terms that change the character of the amplitude significantly and can not be motivated 
physically. 

In accordance with the exact validity of the eikonal formula, the physical situation is 
such, that the transferred momenta are restricted by the condition \k 2 \ <C m 2 . They are, 
however, naturally cut off, if one introduces a form factor to account for the finite extent 
of the nuclei. Thus, large momenta have to be cut off at k± w 1 jr > , where r> is the larger 
value of either the nuclear radius or the Compton wavelength of the scattered particle 
||21|| . In this respect, the electron is an exception, since all other particles have Compton 

5 Note, that for the Schrodinger case the exact validity of the eikonal formula can be proven for a 
certain off shell domain of the momentum transfer for the whole energy plane . 
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wavelengths smaller or comparable to the nuclear size. To present the calculations in a 
uniform manner, we use the form factor of the nucleus to cut off the large momenta. 

The impact parameter dependent cross section for particle production, described in 
the equivalent photon method reads |L9| 



duo 



n|| (wi, uz, b)ar(u 1 ,u 2 ) + n±{u 1 ,u i , 6)<tX 7 (^i, u 2 ) 



(15) 



with the two-photon distribution functions n\\/±(ui,U2,b). The elementary two-photon 
cross sections and the two-photon distribution functions explicitly account for the parallel 
or orthogonal orientation of the photon-polarizations, denoted by the indices || and _L, 
respectively. Since the integration over the impact parameter plane implies an averaging 
over the photon polarizations, the explicit occurrence of the photon polarizations in the 
impact parameter dependent cross section is expected. The functions n\i/±(ui,u 2 ,b) can 
be expressed in terms of single-photon distribution functions n(uj,b), depending on the 
transverse separation: 



with 



n\\(ut,u 2 ,b) 
n±(u u u 2l b) 

n(uj, b) 



d 2 x±n(uji, x± — b) n(ui2, x±) 
d 2 x±n(uji, x± — b) n(uj 2 , x±) 



(x±_ — b) ■ x± 
\x± — b\ \x±\ 

(xj_ — b) x x± 



\x± 



F-L 



(16) 
(17) 



Z 2 a 
n 2 ui 



Ji is a Bessel function. The function F denotes the chosen form factor of the nucleus. 

For a pointlike charge (F = 1), the photon distribution function can be calculated 
analytically. We obtain for the Coulomb potential and the modified potential 



n(u, b) 



Z 2 a u 

7J-2 ry2 



Z 2 a u; 



Ujb" 

, 7 , 



K 



1+iZa 



retarded Coulomb potential 



modified potential (6) 



(19) 



V{l-iZa) 

K v is a modified Bessel function. For small arguments of the Bessel function one can use 



the asymptotic expression [22 



(20) 



Therefore, for ub <C y and assumed point-like charges the photon distribution functions 
(|19|) nearly completely agree. 
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Figure 2: The single-photon distribution function for various photon energies (indicated 
in the plot) as a function of the transverse distance from the ultrarelativistic charge. The 
calculations are done for lead ions (Z = 82) at LHC energies (7 ~ 3000). 

We have numerically evaluated the photon distribution function for an extended 
nucleus, using a gaussian form factor F(Q 2 ) = e - Q ' 2 ^ 2Q ^ with Q = 60MeV which de- 
scribes the Pb nucleus |23| . Figure ^| shows a comparison between the photon distribution 
functions for both, point like nuclei (F = 1) and extended nuclei usin either the retarded 
Coulomb potential or the modified potential (H), which represents the exact calculation. 
At small distances up to a few multiples of the nuclear radius, the modified potential 
gives a smaller number of equivalent photons than the pure Coulomb potential. For large 
distances or large photon energies looses its validity and the modified potential gives 
a larger number of photons than the Coulomb potential (see Figure |3|). 

The significance of these results can be assessed in single-photon induced processes, 
such as the electromagnetic dissociation of nuclei in peripheral heavy-ion collisions. In the 
Weizsacker- Williams approach the dissociation probability at a given impact parameter is 
calculated by integrating the product of the measured photon-nucleus dissociation cross 
section of the target, o^{uS), and the equivalent photon distribution, n P (u,b), of the 
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projectile fll8|) over the photon energy |24|, ^5 



P(6) = / duj np(uj,b) a^(uj) 



As shown in [24] for Pb + Pb collisions at LHC energies, the dissociation probability 
violates unitarity at small impact parameters up to 25 fm. At small impact parameters, 
however, the photon distribution of the modified potential is suppressed relative to that 
of the pure Coulomb potential and yields a reduction of the probabiltiy, whereas at large 
impact parameters, the dissociation probability is enhanced. The reduction visible in 
Figure ||] taken alone is too small to cure the unitarity violation. Possibly this can be 
achieved if in addition multi-phonon excitations |27], [26| are taken into account, which 
also reduce the dissociation probability at small impact parameters. 




Figure 3: The ratio of the equivalent photon numbers of the pure Coulomb potential, 
ric(w,x±), and the photon numbers of the effective potential, n e ff(u,x±). At small 
transverse distances, one finds a deviation of up to 70% for small photon energies. Far 
outside the nucleus, the photon distribution functions are determined by (|T9|). For large 
distances one asymptotically finds a deviation in the order of 40%, independent of the 
photon energy. 

The impact parameter dependent two-photon distribution functions also have to be 
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corrected according to ([TJ), thus correcting the pair production probability ([15]), i.e. @, 
for Coulomb effects. Due to the complicated convolution of the single-photon distribution 
functions in (|1~6D and (0), effects on the two-photon distribution function are not obvious. 
For large distances x± ^> b, however, one directly finds an enhancement of the equivalent 
two-photon numbers. 



5 Summary 

It was shown by either summing the perturbation series or by matching plane waves at 
the delta function potential on the light front [ffl [12], that the eikonal expression for the 
scattering amplitude becomes exact in the ultrarelativistic limit (7 — > 00). This allows to 
neglect the squared momentum transfer k 2 relative to the term 2pi ■ k in the denominator 
of the propagator of the scattered particle. As a consequence the applicability conditions 
of the Weizsacker- Williams method are fulfilled automatically - irrespective of the mass 
of the scattered particle. 

Furthermore, the exact validity of the eikonal formula for ultrarelativistic scattering 
processes confirm the Pomeranchuk theorem, stating that the cross sections for antiparti- 
cle and particle scattering at a given target become identical in the ultrarelativistic limit. 
In analogy to pomeron exchange in hadron physics, on can describe the exact interaction 
as the exchange of an effective photon, according to a modified, effective potential given 
by (Q). The cross section, as a peculiarity of the Coulomb interaction, becomes identical 
to the Mott result. The exact pair production cross section in the field of two ultra- 
relativistic colliding (pointlike) ions also reduces to the second order perturbative result 
12 1 which was evaluated in [14]]. This allows for two conclusions: i) The production rate 



scales with the square of the target and the projectile charge || |10| ii) Asymmetries in 
the electron and positron spectra should not occur. 

Note, however, that the presented formalism is valid only if the produced particles are 
fas with respect to both nuclei. Therefore, the observed || asymmetry at small electron 
and positron momenta remains unaffected by these considerations. 

We applied the Weizsacker- Williams approach to pair production using the modified 
potential @, correctly accounting for the Coulomb boundary conditions. The impact 
parameter dependent single-photon distribution calculated with the modified potential 
shows deviations from the equivalent photon distribution function obtained from the 
retarded Coulomb potential in order of up to 70% at small separations and approximately 
40% at large separations from the ion. 

In combination with multi-phonon excitation, the correct treatment of Coulomb dis- 
tortions possibly solves the problem of unitarity violation in photodissociation processes 
in ultrarelativistic heavy-ion collisions. 

The pair production probability is also subject to changes due to the modified photon 
numbers at given impact parameters and photon energies. The perturbation theoretical 
probability, as calculated here, rather represents the average number of produced pairs 
and exceeds unity at sufficiently small impact parameters. The "true" pair production 
probability has to be corrected by the vacuum-to-vacuum amplitude, which in turn can 



be calculated from the perturbative pair production probability [28, 29]. This nontrivial 



influence on the pair production cross section is subject of further studies. 
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